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ABSTRACT 

The Triangulum-Andromeda stellar clouds (TriAndl and TriAnd2) are a pair 
of concentric ring- or shell-like over-densities at large R (~ 30 kpc) and Z 
-10 kpc) in the Galactic halo that are thought to have been formed from the 
accretion and disruption of a satellite galaxy. This paper critically re-examines 
this formation scenario by comparing the number ratio of RR Lyrae to M giant 
stars associated with the TriAnd clouds with other structures in the Galaxy. The 
current data suggest a stellar population for these over-densities (/rr : mg < 0.38 
at 95% confidence) quite unlike any of the known satellites of the Milky Way 
(/rr:MG ~ 0.5 for the very largest and /rr : mg >> 1 f° r the smaller satellites) and 
more like the population of stars born in the much deeper potential well inhabited 
by the Galactic disk (/rr : mg < 0.01). N-body simulations of a Milky-Way-like 
galaxy perturbed by the impact of a dwarf galaxy demonstrate that, in the right 
circumstances, concentric rings propagating outwards from that Galactic disk can 
plausibly produce similar over-densities. These results provide dramatic support 
for the recent proposal by Xu et al. (2015) that, rather than stars accreted from 
other galaxies, the TriAnd clouds could represent stars kicked-out from our own 
disk. If so, these would be the first populations of disk stars to be found in the 
Galactic halo and a clear signature of the importance of this second formation 
mechanism for stellar halos more generally. Moreover, their existence at the very 
extremities of the disk places strong constraints on the nature of the interaction 
that formed them. 

Subject headings: stars: variables: RR Lyrae - Galaxy: halo - Galaxy: stellar 
content - Galaxy: structure - techniques: spectroscopic - Galaxy: disc. 
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1. Introduction 

Formation scenarios that have been proposed for the diffuse stellar halos that surround 
galaxies can be broadly divided into three categories. First, stars might be accreted from 
other dark matter halos that have been disrupted by the main galaxy (e.g., Searle & Zinn 
1978). The last decade has revealed abundant evidence for this mechanism in the form of 
debris structures in space and velocity around our own (e.g. Majewski et al. 2003; Belokurov 
et al. 2006; Schlaufman et al. 2009) and other (e.g. M31, see Ibata et al. 2007; Gilbert et al. 
2009) galaxies. Second, stars might form in situ from gas located in the the halo itself (e.g. 
Eggen et al. 1962). While such populations have been seen in hydrodynamic simulations 
of galaxy formation (Tissera et al. 2012) it is as yet unclear whether their formation is 
an artifact of either numerical techniques or the assumed star-formation prescriptions. The 
observational evidence for the existence of this population remains controversial (e.g. Carollo 
et al. 2008; Schonrich et al. 2014, for opposing views). Lastly, stars initially made within a 
galaxy’s disk might be kicked-out onto halo orbits by some perturbation (Zolotov et al. 2009, 
2010; McCarthy et al. 2012). Deep star count surveys of M31 have revealed a disturbed disk¬ 
like structure extending to 40-50 kpc from the center of the galaxy (Ferguson et al. 2002) 
which contain stellar populations (Richardson et al. 2008; Bernard et al. 2015) indicating 
that kicked-out stars do indeed contribute to its inner halo (and see Purcell et al. 2010, for 
comparison simulations). However, observations of kicked-out stars around our own Galaxy 
are at this point suggestive rather than conclusive (e.g. Sheffield et al. 2012; Hawkins et al. 
2015) 

The Triangulum-Andromeda clouds are examples of stellar over-densities in the halo 
of the Milky Way thought to have formed through the disruption of a satellite and hence 
contributing to the accreted component of the stellar halo. A single over-density was orig¬ 
inally identified in M giant stars selected from the Two Micron All Sky Survey (2MASS) 
covering the region 100° < l < 160° and —35° < b < —15° at approximately 20 kpc from 
the Sun (Rocha-Pinto et al. 2004). Spectroscopic follow-up of these M giants showed them 
to have a velocity distribution with a small dispersion (^25 km s _1 ) and a shallow gradient 
with mean velocities increasing with decreasing l (in the Galactic-Standard-of-Rest - GSR - 
frame). A double sequence of main-sequence stars in the foreground of the PANDA’s survey 
was later found in the same region (called TriAndl and TriAnd2 by Martin et al. 2007), and 
Martin et al. (2014) subsequently presented evidence for the presence of several additional 
structures. In our own, recent study (Sheffield et al. 2014), we found that the 2MASS data 
suggests that there are M giant counterparts to both of the primary TriAndl and TriAnd2 
sequences. Requiring both M giant and main-sequence stars to belong to the same popula¬ 
tions and fitting isochrones to both simultaneously gave (distance, age, [Fe/H]) estimates of 
(~18 kpc,6-8 Gyrs, -0.75) and (~28 kpc, 10-12 Gyrs, -0.85) respectively (i.e. Galactocentric 
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(R,Z) ~ (24,-8) kpc and (33,-12) kpc for the central pointing at (l,b) = (125°,—25°), 
although the center of the over-density itself could be hidden at lower \z\ behind the Galactic 
disk). Using a spectroscopic survey we showed that M giants in both TriAndl and TriAnd2 
followed the same velocity distribution found in the original study (Figure 1, top panels). We 
proposed that the large Galactocentric radius and distance below the mid-plane combined 
with the small dispersion was suggestive of neither a disk nor random halo population. We 
used simulations of satellite disruption to show that both TriAndl and TriAnd2 could plau¬ 
sibly be due to single accretion event of a satellite on a retrograde orbit — the two distinct 
structures overlapping in orbital phase, but corresponding to material lost on two separate 
pericentric passages. 

In this paper we re-examine our proposal for the origin for TriAndl and TriAnd2 using 
the results from our own spectroscopic survey of RR Lyrae stars in the region. The original 
aim of this survey was to use the spectra to identify RR Lyrae that followed the same velocity 
sequence as the M giants and hence could be associated with these structures. The properties 
of the associated RR Lyrae would provide better constraints on the distances to TriAndl 
and TriAnd2 as well as shed light on their stellar populations. 

In Section 2 we describe our target selection, data collection and reduction. In Section 
3 we present the results of the survey and discuss the implications for the stellar populations 
(as represented by the number ratio of RR Lyrae to M giants stars, /rr ; mg) present in the 
TriAnd clouds. In Section 4 we discuss the implications of our findings the origin of TriAndl 
and TriAnd2 and in Section 5 we present a numerical simulation to illustrate this possible 
formation mechnism. We summarize our conclusions and outline future work in Section 6. 


2. Data 

TriAndl was originally identified as an over-density in M giant stars selected from the 
Two Micron All Sky Survey. Based on the stellar populations of Milky Way satellites and of 
other halo over-densities (e.g., the Hercules-Aquila cloud; Simion et al. 2014), we also expect 
to find an over-density of RR Lyrae stars (RRLs) in this region. RR Lyrae stars have several 
properties that make them useful as tracers of halo over-densities. First, they are bright stars 
(M r = 0.6 mag at [Fe/H] = —1.5 dex) that can be detected at large distances (5-120 kpc for 
14 < r < 21). Second, distances to RRLs measured from optical data are precise to ~ 6% 
(Sesar et al. 2013) versus, e.g., 15-20% for M giants (Sheffield et al. 2014). Finally, RRab 
stars have distinct, saw-tooth shaped light curves which make them easy to identify given 
multi-epoch observations (peak-to-peak amplitudes close to 1 mag in the r-band and periods 
of 0.6 days). Our original goal was to obtain spectra of RRLs selected from the TriAnd 
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region in order to identify stars with velocities consistent with the M giant population, thus 
allowing for an accurate measurement of the distance to the TriAnd structures. 


2.1. Obtaining velocities for RR Lyrae in the TriAnd region 

2.1.1. Target selection: PTF 

We selected RRLs from the Palomar Transient Factory (PTF; Law et al. 2009) in the 
region around TriAndl and TriAnd2: Galactic longitude and latitude, ( l,b ), in the ranges 
100° < l < 160°, —35° < b < —15°, and heliocentric distances in the range 15 < d < 35 kpc. 
We follow the procedure of Drake et al. (2013) and use the absolute-magnitude-metallicity 
relation of Catelan & Cortes (2008) to compute the distances, where we assume all stars 
have a metallicity of [Fe/H] = —1.5. There are 70 RR Lyrae detected in PTF that fall in 
this volume. 


2.1.2. Data collection and reduction 

We used the 2.4 meter Hiltner telescope at MDM Observatory on two observing runs 
(August 2013 and October 2013) to obtain spectra for 20 of the RR Lyrae stars selected from 
PTF. We used the Modular Spectrograph (ModSpec) with the 600 1 mm -1 grating blazed at 
4752 A, centered near the prominent Hct absorption feature, and the Templeton CCD. This 
set-up provides spectral coverage from 5300 A to 7200 A; this range was selected to focus 
on the Ha line, which we use to derive the systemic velocities of the RR Lyrae stars. We 
found that total exposure times of 1200 seconds are adequate to achieve a signal-to-noise 
of S/N « 15 near the Ha line for an RR Lyrae with a V -band magnitude of 17. HgNe 
lamps were taken before and after each observation, and 2 sets of 10 flat field frames were 
taken, one set at the beginning of the night and the other set at the end. Heliocentric radial 
velocities for the Ha line were found using standard IRAF tasks. The motion of the observer 
was corrected for using the rvcorr task. To check for systematic shifts in the wavelength 
calibration, we used the [O 1] night sky emission line at 6300.304 A; no systematic shifts 
were found, although we did shift any spectra that had a deviation greater than 0.1 A so 
that the observed value of the line matched the lab value. 

To derive the systemic velocities of the RR Lyrae stars, we use the procedure outlined in 
Sesar (2012). The systemic velocities are found using the H Balmer lines, although a different 
template is used for each line — we use the template derived for the Ha line. For most stars, 
we took spectra at 2 phase points within the phase range 0.05-0.85, however several stars 
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have only 1 phase point measured. We use standard least-squares optimization to fit the Ha: 
template to the individual radial velocity measurements and derive the systemic velocity for 
each star. We compute the uncertainty in the systemic velocity using Equation 9 of Sesar 
(2012), which assumes all distributions are Gaussian but accounts for the uncertainties in the 
individual velocity measurements and inherent uncertainties in the RR Lyrae magnitude-RV 
amplitude relation and template parameters. We take the systemic velocity to be the radial 
velocity of the fit template at phase 0.27, the canonical phase corresponding to zero pulsation 
velocity. The systemic velocity uncertainties are typically around 15-20 km s . Figure 2 
shows example template fits for two of the observed RR Lyrae stars (top panels) along with 
phase-folded, R-band light curves from PTF (bottom panels). 


2.2. Existing Catalina Sky Survey observations 

The TriAnd region defined in Section 2.1 is also covered by the Catalina Sky Survey 
(CSS; Larson et ah 2003). One additional star in this area has a radial velocity measured by 
CSS follow-up efforts (Drake et ah 2013) and we include this star in our sample of TriAnd 
RR Lyrae velocities. 90% of the PTF RR Lyrae appear in the CSS sample, but there are 
many more CSS RR Lyrae in total. The CSS fields in this region have more observations 
and thus the completeness of RR Lyrae detection of the CSS is better than that of PTF. 
There are 142 RR Lyrae in CSS in the TriAnd volume, compared to 90 from PTF. 


3. Results I: the density and velocity distribution of TriAnd 

We expect to find over-densities of RRLs in the TriAnd region consistent with the M 
giant population. In this section we compare the observed and the predicted number of RRLs 
and evaluate whether an over-density of RRLs is present. We then model the velocities of 
the observed RRLs and find that the RR Lyrae velocities appear to be consistent with a 
pure background (halo) velocity distribution. 


3.1. Density of RR Lyrae stars 

We estimate that the average CSS field in the TriAnd region has ~300 observations. 
This is comparable to the average number of CSS observations along the SDSS Stripe 82 
region, for which a previous study has identified nearly 100% of the RR Lyrae stars in this 
region (Sesar et al. 2010). To estimate the completeness of RR Lyrae identification in CSS 
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as a function of distance, we compare the number RR Lyrae in CSS relative to the number in 
the Stripe 82 sample in several distance bins ranging from 5 to 35 kpc. We find that between 
15-21 kpc, CSS is 80% complete for fields with several hundred observations; there are 77 
RR Lyrae in this distance range, or, a completeness-corrected number TVcss = 77/0.8 ~ 96. 

We compute the predicted number of RR Lyrae stars using the smooth halo density 
relation of Sesar et al. (2011), 

prrO) = po (!) 

p 0 = 5.6 kpc~ 3 (2) 

a = 2.42 (3) 

where r is the Galactocentric radius (with halo oblateness, q = 0.63, taken into account 

Sesar et al. 2011) and we fix the distance of the Sun to the Galactic center, doc — 8 kpc. We 
find the expected number of RR Lyrae stars to be N exp & 225, over twice the completeness- 
corrected number of RRLs measured by the CSS. This is not surprising given the vast amount 
of substructure observed in halo RR Lyrae density (e.g., Sesar et al. 2010, 2013), but this 
does indicate that the region around TriAndl is under-dense , not over-dense, in RRLs. 

We now also compare the smooth model predictions for nearby regions of the sky with 
the number counts of RRLs from CSS. We find that in the adjacent region towards the 
Galactic anticenter, (160° < l < 220°, —35° < b < —15°), the single power-law model 
again vastly over-predicts the number of RRLs at distances between 15-21 kpc; we find a 
completeness corrected number of RR Lyrae in this volume ^55 compared to the predicted 
number ^190. In the opposite direction, the region adjacent and closer to the Galactic center 
(but slightly lower in b due to survey footprint, 40° < l < 160°, —40° < b < —20°) contains 
a completeness corrected estimate of ^350 RRLs in the CSS and a predicted number ^360. 
Again, in these high (absolute) latitude fields, many of the stars are halo stars where other 
authors have seen a large amount of substructure and variance to the RR Lyrae number 
counts (e.g., Sesar et al. 2010, 2013). 


3.2. Analysis of the line-of-sight velocity distributions 

To quantify the fraction of RR Lyrae in our sample that may belong to the TriAndl 
feature, we develop a statistical model and simultaneously fit the velocity profile and disper¬ 
sion of the TriAnd M giants from Sheffield et al. (2014) with the RR Lyrae stars to assess 
how many RR Lyrae may be associated with the cold sequence. 
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We model the line-of-sight velocities using a Gaussian mixture model: any individual ve¬ 
locity lias some probability of being drawn from either an isotropic halo velocity distribution, 
or from a narrower TriAnd sequence: 

p(v | $ TriAnd; /) = (!- f)N(v\v h ,crl) + fp(v I ^TriAnd ) (4) 

where Uh and <Th are the mean halo velocity and the halo velocity dispersion — fixed to 
(uh, cr h ) = (0,106) km s' 1 (Brown et ah 2010) — / is the fraction of stars of a given type that 
belong to TriAnd, ^TriAnd are the parameters for describing the TriAnd velocity sequence, 
and J\f(x | /i, a 2 ) is the normal distribution over x with mean p and variance cr 2 . If we assume 
that the TriAnd stars are also Gaussian distributed around a linear velocity trend in Galactic 
longitude, l , we may write 

p(v | ^ TriAnd) = M {v \ V 0 + (dv/dl) l, <7^) (5) 

B TriAnd = (v 0 ,dv/dl,crl) (6) 

where v 0 and dv/dl are the zero-point and velocity gradient of the TriAnd feature, a v is the 
intrinsic velocity dispersion of the TriAnd stars, and we assume the l are measured perfectly. 

Of course, the line-of-sight velocities are measured with finite uncertainties. If we assume 
the measured GSR velocity, w, is also Gaussian distributed with known uncertainty a w , we 
may write the joint likelihood over measured and true velocity as 

p(w, V | 0 TriAnd; /) = N(w \ V, a^) p(v | ^TriAnd, /)• (7) 

Marginalizing over the true velocity, v, we are left with the per-star likelihood of measured, 
GSR velocity, w, in terms of the mixture model parameters: 

p{w | $ TriAnd; /) = J p(w, V \ ^TriAnd, /) dv ( 8 ) 

= (1 - /) A f(w I Uh, <?l + cTh) + fAf(w I v 0 + (dv/dl) l, a 2 w + a 2 ). (9) 

The M giant and RR Lyrae samples (with observed GSR velocities wmg and wrrl, respec¬ 
tively) will have different TriAndl membership fractions, /mg and /rrl- We may therefore 
write the full likelihood for both populations as 

p(v | 6 TriAnd; /mG; /rrl) = p(v MG | ^TriAnd; /mg) P(wRRL | ^TriAnd; /rrl) (10) 

where the underlying model for the velocity sequence, parametrized by ^TriAnd; is the same. 
Assuming the velocity measurements are independent, the likelihood over all measured ve¬ 
locities, {wi}, is the product of the individual (marginal) likelihoods: 

C = p({Wi} | 0TriAnd; /mG; /rrl) = ]^[ p{w l \ ^TriAnd; /mG; /rrl)- (H) 


The posterior probability is then 

p{&' TriAnd, /mG, /rRL | {Wi}) = ^ P{{Wi} | ^TriAnd, fuG, /rRl) ^(^TriAnd) p(/mg) P(/rRl) (12) 

where the evidence integral, Z, is a constant that depends only on the data, p(^TriAnd) is the 
prior probability over the TriAnd velocity model parameters, and p(f mg) and p(/rrl) are 
the prior probabilities over the fraction of M giants and RRLs, respectively, that are halo 
(field) stars. We assume a uniform prior for the zero-point of the TriAnd velocity sequence, 
Vo, over the range (—200,200) km s -1 , a scale-invariant (logarithmic) prior over the velocity 
dispersion of TriAnd, <7%, and assume uniform priors over each membership fraction over the 
range (0,1). 

We use an ensemble Markov Chain Monte Carlo (MCMC) algorithm (Goodman & 
Weare 2010) implemented in Python (Foreman-Mackey et al. 2013) to draw samples from 
the posterior probability distribution defined above (Equation 12). The algorithm uses an 
ensemble of “walkers” that naturally adapt to the geometry of the parameter-space; we use 
32 walkers in this analysis. We burn-in the walkers for 256 steps starting from randomly 
drawn initial conditions (sampled from the priors) and then re-initialize the walkers from 
their final burn-in positions and run for another 4096 steps. Figure 3 shows contour-plot 
projections of the samples from the posterior. We find that the TriAnd velocity sequence 
has a velocity dispersion, gradient, and zero-point of 


<j v = 22.0 ± 3.6 km s 1 


(13) 

(14) 

(15) 


dv/dl = —0.86 ± 0.24 km s 1 deg 1 
Vo = 155 ± 32 km s -1 


but with significant correlation between Vo and dv/dl. 

For each star (M giant and RR Lyrae), we compute the the posterior probability that 
the star is a member of the TriAnd velocity sequence (Foreman-Mackey 2014). To do this, 
we make use of an additional binary, nuisance parameter, q — when q = 0, the star is part 
of the halo population, and when q — 1, the star is a member of TriAnd. The conditional 
probability that a given star is a member of TriAnd (e.g., q — 1) is 


p(q= 1 | ^TriAnd, f,w) 


f N{w I Vo + (dv/dl) l, aj + a 2 v ) 
p(w | # TriAnd, /) 


(16) 


and the posterior probability is just the integral of this expression over all parameters. 
This integral can be approximated with a Monte Carlo integral using the posterior samples 
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generated with MCMC, 

p(q = 1 | w) = j p(q — 1 | 0TriAnd, /, w) p{0 T riAn d, f\w)df d0TriAnd (17) 

1 N 

TWnd = ! I 0%L d , f (n) ,w)- (18) 

n 

Figure 1 (top panel) shows the M giant stars with -PrriAnd > 0.8 in black, and M giants 
identified as being background halo stars in grey. 


3.3. Estimated Ratio of RR Lyraes to M giants in TriAnd: /rr : mg 

We simultaneously constrain the fraction of RR Lyrae that may belong to the cold 
TriAnd velocity sequence with our velocity model. Figure 3 (bottom right) shows the 
marginal posterior distribution over /rrl; the MAP value of /rrl = 0.004, and the val¬ 
ues /rrl = (0.183,0.38) contain 68% and 95% of the probability mass, respectively. We use 
the marginal posterior over /rrl, the total (completeness-corrected) number of RR Lyrae 
in this region, Acss, and the number of M giants identified by Sheffield et ah (2014) to be 
in TriAndl, IVmg — 74, to compute the distribution over the ratio of the number of RR 
Lyrae to the number of M giants in TriAnd, /rr : mg = /rrl ATcss/^Yivig- Figure 4 shows our 
posterior constraints on this ratio. 

To validate our model and test whether we would expect to detect a nonzero fraction 
of RR Lyrae stars associated TriAnd with our data, we generate mock RR Lyrae velocity 
data and use the model to infer /rr : mg f° r many trials. We assume that there are 100 RR 
Lyrae in total and 75 M giants in our mock TriAnd region, but only “observe” 20 RR Lyrae 
stars assuming 40% of the stars are sampled from a low-dispersion, linear velocity sequence 
with an intrinsic scatter set to 20 km s -1 . These numbers produce a true RR Lyrae to M 
giant fraction, /rr : mg, equal to ~0.5, which represents a lower bound for what we would 
expect for accreted stellar populations (see Section 4). The rest of the stars are drawn from 
the same isotropic halo distribution used above. For each mock RR Lyrae star, we observe 
the velocity with an uncertainty drawn from a uniform distribution over the range 15-20 
km s -1 . For this test, we assume we know the true TriAnd velocity trend and only consider 
a single free parameter: the fraction of RR Lyrae that belong to the cold sequence, which we 
convert to /rr-mg using the numbers above. We perform 1024 trials in this test. Figure 5 
shows histograms of posterior samples over /rr : mg f° r each of the 1024 trials, plotted with 
transparency. The green line indicates the true fraction in the mock dataset. The posterior 
distribution over /rr : mg f° r any individual trial with a true ratio of /rr : mg ~ 0.5 generally 
recovers a fraction between 0.3-0.7 and is almost never peaked towards zero. 
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4. Results II: Implications for the origin of the TriAndl and TriAnd2 

substructures 

4.1. Accreted satellite? 

The implications of not being able to easily find RR Lyrae associated with TriAnd conlcl 
be very interesting. All dwarf galaxies that are close enough to have been surveyed for RR 
Lyrae are found to contain at least a few (Clementini 2010). However, only the largest 
satellites (e.g. the Small and Large Magellanic Clouds - SMC, LMC - and the Sagittarius 
dwarf spheroidal galaxy - Sgr) are known to contain M giants. This systematic difference 
between low- and high-luminosity dwarf galaxies is easily explained as a consequence of the 
age and metallity distributions of their stellar populations. Single age populations old enough 
to containing stable, helium core-burning stars have to be less metal rich than [Fe/H] ~ —1 
for their horizontal branches to extend beyond the instability strip that RR Lyrae stars 
inhabit. Conversely, populations have to be more metal rich than [Fe/H] « —1 for red- 
giant-stars to evolve through spectral types as late as M. Only the LMC, SMC, and Sgr 
have metallicity distributions that are both metal-rich enough to contain M giants yet with 
a sufficiently populated tail to lower metallicities that they also contain significant numbers 
of RR Lyrae (e.g. Chou et al. 2007; Carrera et al. 2008; Kirby et al. 2011). (Note, while the 
ages and metallicity estimates from isochrone fitting in Sheffield et al. 2014, are inconsistent 
with RR-Lyrae populations, these estimates were made under the assumption that the main- 
sequence detections corresponded to the overdensities seen in the M giant populations and 
hence do not themselves rule out the presence of RR Lyrae in these structures.) 

If the TriAnd structures did form from the disruption of a satellite, then the presence of 
significant M giants in them suggests that the progenitor was likely to have stellar populations 
like the LMC, SMC, and Sgr. Hence, we can use the ratio of the numbers of RR Lyrae and 
M giants, /rr : mg, in these satellites as a guide for how many RR Lyrae we might find also 
associated with these structures if they were likely accreted from a large satellite. 

Figure 6 summarizes how we estimated /rr : mg f° r Sgr. The top left-hand panel shows 
the positions of RR Lyrae in the SAG plate of the catalogue of Cseresnjes et al. (2000), 
which maps the outer structure of Sgr and its debris at Galactic latitudes above the main 
progenitor, —14° < l < —5°. Only those with R-magnitudes greater than 18.5 are shown, 
which selects Sgr stars with only mild contamination from the bulge (see Cseresnjes et al. 
2000, Figure 3). The top right-hand panel shows the M giants in same region selected from 
the 2MASS catalogue (following Majewski et al. 2003) in the same areas. The M giants were 
chosen within appropriate color limits (after reddening correction and to match the selection 
criteria, 0.9 < (J — K s ) < 1.3, employed for the TriAnd sample in Sheffield et al. 2014) 
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and chosing (extinction-correction) magnitudes (8 < K s < 13) that would encompass the 
full range possible given the distance to Sgr and the known spread in M-giant luminosities. 
The bottom left-hand panel repeats the middle panel, but with M giants with reddening- 
estimates E (B — V) > 0.2 highlighted in red. Clearly, the red points in this panel cluster 
around the apparently “empty” regions of the plot, and the clumpy nature of the M giant 
distribution compared to the RR Lyrae can be attributed to reddening effects. To address this 
incompleteness we placed a grid of overlapping fields of radius 0.5° over the area and made 
separate estimates of the number of M giants and RR-Lyrae for only those that contained 
no M giant with E (B — V) >0.2 (highlighted by blue points in the same panel). 

The bottom right-panel of Figure 6, which shows the magnitude-distributions for M 
giants for the fields highlighted in blue in the bottom left-hand panel, illustrates how the 
final estimates for /rr : mg were made. For the M giants, the broad spread in absolute 
magnitudes causes significant overlap in the apparent magnitude distributions of the (brigter, 
K s ~ 9 — 10) bulge stars and (fainter, K s ~ 11 —12 ) Sgr stars, so in this case the contribution 
of Sgr stars in a given field was estimated by fitting two Gaussians (with means, dispersions, 
and amplitudes as free parameters) to the full distribution and estimating the number of 
Sgr stars in the combined fields from the ratio of the the amplitudes. Overall, we found a 
total of 2319 RR Lyrae and 4770 M giants across all fields associated with Sgr, suggesting 
/rr:MG ~ 0.5. The ratio in the 41 (overlapping) individual fields had a mean and dispersion 
i n /rr:MG °f 0.48 and 0.088 (and varied across the full range 0.31-0.69). The individual 
fields had ~ 50-100 and ~ 100-200 in them respectively, suggesting that fluctuations in 
counts should lead to a dispersion in /rr : mg of order 0.02. We interpret the larger measured 
dispersion as due to a combination of possible factors including intrinsic variations in stellar 
populations, uncertainties in the counts themselves in these crowded fields and Gaussians 
not being an appropriate form for the fitting function. 

Figure 7 repeats the same plots as Figure 6 for the LMC, where the RR Lyrae counts 
come from the catalogue of the OGLE team (Soszynski et al. 2009). The same approach as 
for Sgr was used to find an overall number of 2686 RR Lyrae and 5210 M giants in the fields 
we examined, suggesting /rr : mg ~ 0.5. The ratio in the 18 (overlapping) individual fields 
had a mean and dispersion in /rr : mg of 0.52 and 0.086 (and varied across the full range 
0.39-0.69). The individual fields had ~ 100-200 RR Lyrae and ~ 200-400 M giants in them. 
We again interpret the large dispersion as due to a combination of the same possible factors. 

The estimates and ranges for /rr : mg from our analyses of Sgr and the LMC are plotted 
as lines and shaded regions on top of the PDF for /rr : mg generated from our data in the 
TriAnd region in Figure 4. While an Sgr- or LMC-like population cannot yet be absolutely 
ruled out, the data suggests that the TriAnd structures (/rr : mg < 0.38 at 95% confidence) 
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have a much lower ratio of RR Lyrae to M giants than these objects (/rr : mg ~ 0.5 ±0.1), 
and very unlike any of the other known satellites which are not known to contain M giants 
at all. 

Of course, the disrupted object might not have populations like the Sgr or LMC. However 
a lower luminosity or more ancient satellite would only have older and/or lower metallicity 
populations, suggesting a larger fraction of RR Lyrae. A more massive satellite might be 
higher metallicity overall and contain a smaller fraction of RR Lyrae, but then would be 
likely to give TriAnd a larger dispersion than observed and leave many other signatures of 
its recent disruption in the disk and the halo. Overall, we conclude that the lack of RR 
Lyrae stars contributing the cold velocity sequence outlined by M giants in the region calls 
into question the interpretation of the origin of TriAnd structures as satellite debris. 


4.2. Kicked-out disk? 

In comparison to the Galactic satellites, the Galactic disk is more metal rich and its 
metal-poor tail is much more poorly populated (Schlesinger et ah 2012) — hence it contains 
a much higher fraction of M giants. The observed number density for RR Lyrae at the 
disk midplane around the Sun is ps 5-10 per kpc 3 (predominantly from the thick disk, see 
Layden 1995; Amrose & Mckay 2001) while the number density of M giants predicted by 
the Galaxia stellar populations code is ~ 4000 per kpc 3 (Sharma et al. 2011), suggesting 
a number ratio of less than 1%. Hence, the apparent lack of RR Lyrae associated with 
the TriAnd velocity sequences could be indicative of these populations originating in our 
Galactic disk. The ages and metallicities of the M giants in these sequences indicated by 
isochrone fitting in Sheffield et al. (2014) (ages of 6-8Gyrs vs 10-12 Gyrs, and [Fe/H] ~ -0.75 
vs - 0.85 for TriAndl and TriAnd2 respectively) are loosely consistent with an extrapolation 
of the scales and gradients of metallicities seen at smaller radii (3-15kpc) in the disk, along 
with the models of age distributions that these observations suggest (Hayden et al. 2015). 
However, there is not yet enough known about stellar populations at these extreme radii for 
this extrapolation to be conclusive in itself. 

We conclude that the stellar populations of TriAndl and TriAnd2 suggest these sub¬ 
structures are likely due to a perturbation kicking stars out of the disk rather than a satellite 
accretion event. This formation scenario has not been considered until very recently (Xu 
et al. 2015) due to the large Galactocentric radius (30-40 kpc) and position below the disk 
plane (5-10 kpc) of these features. In the next section we examine whether coherent stellar 
rings can be kicked to such extreme locations through encounters with known satellites. 
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5. Results III: the creation of extreme disk populations 

In order to investigate the plausibility of creating TriAndl and TriAnd2 with a dynam¬ 
ical disturbance to the Galactic disk, we performed N-body simulations of three-component 
(disk, bulge and dark-matter halo) Milky-Way-like galaxies being perturbed by a satellite 
with the mass and orbital characteristics of the Sagittarius dwarf galaxy. The initial con¬ 
ditions were generated using the rejection method (Kuijken & Dubinski 1994) to sample 
the distribution functions of each components as in Widrow et al. (2008). This method is 
reliable and shown to produce more stable galaxies (Kazantzidis et al. 2004) than methods 
relying on moments of the collisionless Boltzmann equation which have typically been used 
in the literature (Hernquist 1993). The simulations were evolved with the Gadget-3 treecode 
(Springel et al. 2001; Springel 2005). 

There exists prior work in the literature on the impact of MW satellites (LMC or Sgr) 
on our Galactic disk(e.g. Weinberg 1995; Bailin 2003). In particular, (Purcell et al. 2011) 
qualitatively show that the impact of a dark matter halo of the size that could have hosted 
the progenitor of Sgr could plausibly produce the Monoceros Ring or the Giant Anticenter 
Stellar Structure (Yanny et al. 2003; Crane et al. 2003) as well as the north/south velocity 
asymmetries in the Galactic disk that have been seen in the RAVE and SEGUE surveys 
(Widrow et al. 2012; Gomez et al. 2013; Williams et al. 2013). As a starting point to our 
investigation, we have re-run a similar N-body simulation to Purcell et al. (2011). The dark 
matter halo of the MW galaxy is given by an NFW halo with a virial mass of M V i r = 10 12 M 0 
and a scale radius of r s = 14 kpc. The disk has a mass of 3.7 x 10 10 M 0 , an exponential 
scale length of 2.8 kpc and a vertical scale height of 0.4 kpc and a cut-off radius at ~ 50 
kpc. The bulge has a mass of 9.4 x 10 9 M 0 and a Sersic profile with index n = 1.28 and 
effective radius of 0.56 kpc. The Sagittarius analogue is modeled with a virial mass of 
M vir = 10 11 M 0 and scale length r s = 6.5 kpc. The satellite is launched at 80 kpc from 
the Galactic centre from the plane of the MW on a vertical trajectory with a velocity vector 
v = (v x ,v y ,v z ) = (0, 0, —80) km s -1 towards the Galactic North. At those distances, the 
halo of the satellite galaxy will have experienced significant mass-loss due to stripping so we 
artificially take into account this effect by truncating the density profile at the tidal radius 
rt = 30 kpc. This makes the satellite mass a factor of m3 smaller than the initial virial 
mass originally set. The simulation is run for 2.4 Gyr and the Sgr dwarf reaches present-day 
configuration at m 2.3 Gyr. Our simulation followed the evolution of m 7 x 10' particles 
(-Wisk = 6 x 10 6 , iVbuige = 2 x 10 6 , A r haio = 4 x 10', iVgg r = 4.2 x 10 6 ). The dark matter particle 
mass was ~ 4 times larger than disk particles. The initial disk model is tested for stability by 
running it for a period of 2 Gyr. Apart from the appearances of transient flocculcnt spirals 
the disk showed no variations in its density profile and showed no appreciable disk thickening 
due to interactions with the more massive MW surrounding dark matter particles. 
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Figure 8 illustrates the results at the end point of 2.3 Gyr of evolution in our simulation. 
The top-left, bottom-left, and bottom-right panels show two-dimensional maps of the over¬ 
density, average z-position, and z-dispersion for particles in the disk component. Top-right 
shows an edge-on view of the simulated disk. Outward-propagating ring-like features are 
a natural product of the interaction, with several rings produced during the course of the 
simulation, over-dense by a factor of zs 2 relative to the rest of the disk and spaced by several 
kpc. 

The bottom-right panel demonstrates that these over-densities oscillate in z above and 
below the plane, with the amplitude of the oscillation increasing with R to as much m 4 kpc 
at R = 30 kpc. The thickness of the disk — as indicated by the dispersion ps 0.5-1.0 kpc 
in the top-right and bottom-left panel — also varies. An observer at 8 kpc along the y -axis 
in our simulation would clearly see two of these oscillations as rings several kpc below the 
galactic disk at roughly 18 and 40 kpc beyond their radius. (A third, thinner ring appears at 
~ 10 kpc distance, but its significance is less apparent.) The properties of these simulated 
structures are reminiscent of the GASS and TriAnd clouds (albeit with smaller offsets from 
the disk plane) and support the hypothesis that these structures could plausibly be caused 
by the same interaction. 

There are many things that could explain why the observed oscillations are larger than 
those seen in the simulations. The initial conditions in our simulation — with the satellite 
artificially truncated at its pericentric tidal radius and the disk in equilibrium — neglects 
the first tidal interactions between the two. Moreover, Sgr is not the only object in orbit 
around the Milky Way capable of significantly perturbing the disk: the effect of the LMC 
should also be taken into account. Further, the disk itself is already known to be warped; 
whether this warp is a signature of past or present interactions, its presence will affect any 
further reactions to ongoing perturbations. Finally, the properties of neither the stellar disk 
nor the Galactic potential are well-known in this region and could significantly affect the 
response to any interaction. 

The above list of possible explanations is challenging to investigate. At the same time, 
the presence of this disk population at such extreme Galactic radius and distance from the 
Galactic plane are challenging to reproduce, which suggests that their existence could provide 
strong tests of the various hypotheses. 
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6. Conclusions and Future Work 

By looking at /rr : mg — the ratio of the number of RR Lyrae to M giant stars — in the 
TriAndl structure, we have shown that the stellar populations in this M giant over-density 
in the halo appear more like those in the disk of our Milky Way than in any of the Galactic 
satellites. We have also confirmed, using N-body simulations, that the impact of a large 
satellite can plausibly generate concentric rings of stars in the outskirts of the Galactic disk 
propagating outwards from the Galactic center (albeit with lower amplitude). We conclude 
that TriAndl and TriAnd2 could plausibly be two such propagating stellar rings. As such, 
they would represent the first disk populations seen at such large R and Z around the Milky 
Way. Their existence provides evidence for stars currently being kicked-out of the disk, and 
strong support for this formation mechanism contributing a significant fraction of our stellar 
halo. 


This work bolsters the tentative connections proposed by Xu et ah (2015) between vari¬ 
ous over-densities seen in the star counts in the outer disk and the local velocity asymmetries 
in the disk, which they propose to represent the collective response to an external pertur¬ 
bation. Further work on abundance pattern (as started in Chou et al. 2011) and stellar 
populations could help confirm or rule out this scenario. A similar test of number ratios 
could be used to further examine whether other over-densities at smaller R and Z (such as 
GASS - the Giant Anticenter Stellar Stream - also known as the Monoceros Ring; Rocha- 
Pinto et al. 2003; Ibata et al. 2003) could be the inner rings excited by the same disturbance 
to the Galactic disk. 

Our results also motivate more careful and complete numerical explorations of dynamical 
interactions in order to understand under what circumstances disk material can be kicked 
coherently to such extreme locations. Analogous work on the signatures in outskirts of the 
gaseous Galactic disk (Chakrabarti et al. 2011) has recently led to the identification of a new 
candidate dwarf satellite at 90 kpc from the Milky Way (Chakrabarti et al. 2015). Hence we 
anticipate that the mere existence of the TriAnd structures could place strong constraints 
on such interactions and the potential of the Milky Way. 
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Fig. 1.— Galactocentric radial velocities as a function of Galactic longitude, l, for M giants 
and RR Lyraes in the TriAnd region. Top left panel shows the velocities of all M giants in 
the TriAnd region: black circles are stars with >80% probability of belonging to TriAnd, 
grey circles have <80% probability of belonging to TriAnd and are probably typical halo 
stars. Velocity uncertainties are typically around the size of the marker. Top right panel 
shows a ID histogram of the M giant velocities. Bottom left panel shows a summary of the 
MAP fit of our velocity model (Section 3.2) for the TriAnd sequence (solid grey line), plus 
or minus the derived MAP velocity dispersion of TriAnd (dashed lines). Over-plotted are 
the velocities of our sample of RR Lyrae stars in the TriAnd region. Bottom right shows the 
same M giant histogram (grey), with a histogram of the RR Lyrae velocities (black) with 
the same bins. 
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Fig. 2.— Examples of template fits to radial velocity measurements of two of our sample of 
RR Lyrae stars (top). Black lines show the best-fit template. Grey, shaded region show the 
68% confidence region. Bottom panels show phased, R-band light curves for the two stars, 
demonstrating the data quality of PTF and two different extremes of number of observations. 
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Fig. 3.— Projections of the TriAnd velocity model posterior probability distribution (Equa¬ 
tion 12). Contours show (approximate) 68% and 95% confidence regions in each projection. 
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/RR:MG 

Fig. 4.— Binned samples from the marginal posterior probability over the number ratio of 
RR Lyrae to M giants in TriAnd (histogram). Dashed, vertical line shows the value of the 
ratio typical to the Milky Way disk (<1%). Solid, vertical line and shaded rectangle indicate 
the approximate mean and dispersion of /rr : mg for counts around both Sgr and the LMC 
(see text for exact numbers). 
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Fig. 5.— Histograms of samples from 1024 individual posterior probability distributions 
(plotted with transparency) over the number ratio of RR Lyrae to M giants in mock data 
sets of RR Lyrae velocities. The RR Lyrae velocities are drawn from either an isotropic halo 
distribution or a thin, linear velocity sequence with known properties and intrinsic velocity 
dispersion set to 20 km s -1 . Each velocity is observed with an uncertainty drawn from a 
uniform distribution over the range 15-20 km s -1 . Vertical (green) line shows the true, input 
value of the ratio. 
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Fig. 6.— Top left panel shows the position of RR Lyrae in the SAG plate of the Cseresnjes 
et al. (2000) catalogue. Top right panel shows M giants selected from 2MASS in the same 
region (see text for details). Bottom left panel repeats the top right panel with stars color 
coded in red that have E (B — V) > 0.2 and blue to show fields where stars were counted 
to calculate /rr : mg- The bottom right panel is a histogram of the apparent magnitude 
distribution of the M giants in all the blue fields combined. 
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Fig. 7.— Top left panel shows the position of RR Lyrae in LMC region in the Soszynski et al. 
(2009) catalogue. Top right panel shows M giants selected from 2MASS in the same region 
(see text for details). Bottom left panel repeats the top right panel with stars color coded in 
red that have E (B — V) > 0.2 and blue to show fields where stars were counted to calculate 
/rr:MG- The bottom right panel is a histogram of the apparent magnitude distribution of 
the M giants in all the blue fields combined. 
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Fig. 8.— Top left-hand panel shows the density contrast map of our simulated galactic disk 
at 2.3 Gyr (S(x, y)Z'E(R), where £(i?) is the azimuthally averaged density and £(ic,?/) is the 
local surface density evaluated at position (x,y)). Top right-hand panel shows an edge-on 
view of the disk — vertical oscillations are most apparent towards negative y. The bottom 
panels show the mean position (right) and dispersion (left) of particles perpendicular to the 
galactic disk. The diamond in the bottom right panel indicates the instantaneous position of 
the perturbing satellite (projected) and the lines indicate position of the Sun and the range 
of Galactic longitude occupied by the TriAnd Clouds. 

















